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Abstract. Using recent measurements of angular size of high-z milliarcsecond compact 
radio sources compiled by Gurvits, Kellermann & Frey (1999) and X-ray gas mass fraction 
of galaxy clusters published by Allen et al. (2002,2003), we explore their bounds on the 
equation of state, oi, = PxIPx, of the dark energy, whose existence has been congruously 
suggested by various cosmological observations. We relaxe the usual constraint oj^ > -1, 
and find that combining the two databases yields a nontrivial lower bound on oi,-. Under the 
assumption of a flat universe, we obtain a bound -2.22 < oj^ < -0.62 at 95.4% confidence 
level. The 95.4% confidence bound goes to -1 < ai, < -0.60 when the constraint oi, > -1 
is imposed. 
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1. Introduction 

One of the most remarkable cosmological findings of recent years is, in additional to the cold 
dark matter (CDM), the existence of a component of dark energy (DE) with negative pressure 
in our universe. It is motivated to explain the acceleration of the universe discovered by distant 
type la supernova (SNela) observations (Perlmutter et al. 1998, 1999; Riess et al. 1998, 2001), 
and to offset the deficiency of a flat universe, favoured by the measurements of the anisotropy 
of CMB (de Bernardis et al. 2000; Balbi et al. 2000, Durrer et al. 2003; Bennett et al. 2003; 
Spergel et al. 2003), but with a subcritical matter density parameter ~ 0.3, obtained from 
Send offprint requests to: Zong-Hong Zhu 
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dynamical estimates or X-ray and gravitational lensing observations of clusters of galaxies(for 
a recent summary, see Turner 2002). While a cosmological constant with = -pA is the sim- 
plest candidate for DE, it suffers from the difficulties in understanding of the observed value 
in the framework of modern quantum field theory (Weinberg 1989; CaiToll et al. 1992) and the 
"coincidence problem", the issue of explaining the initial conditions necessary to yield the near- 
coincidence of the densities of matter and the cosmological constant component today. In this 
case, quintessence (a dynamical form of DE with generally negative pressure) has been invoked 
(Ratra and Peebles 1988; Wetterich 1988; Caldwell, Dave and Steinhardt 1998; Zlatev, Wang 
and Steinhardt 1998). One of the important characteristics of quintessence models is that their 
equation of state, = Px/px, vary with cosmic time whilst a cosmological constant remains 
a constant a>A = - 1 • Determination of values of cOx and its possible cosmic evolution plays a 
central role to distinguish various DE models. Such a challenging has triggered off a wave of 
interest aiming to constrain a>x using various cosmological databases, such as SNela (Garnavich 
et al. 1998; Tonry et al. 2003; Harris et al. 2003; Knop et al. 2003; Zhu and Fujimoto 2003); old 
high redshift objects (Lima and Alcaniz 2000a); angular size of compact radio sources (Lima 
and Alcaniz 2002); gravitational lensing (Chae et al. 2002; Sereno 2002; Dev, Jain and Mahajan 
2003; Huterer and Ma 2003); SNela plus Large Scale Structure (LSS) (Perlmutter, Turner & 
White 1999); SNela plus gravitational lensing (Waga and miceli 1999); SNela plus X-ray galaxy 
clusters (Schuecker et al. 2003); CMB plus SNela (Efstathiou 1999; Bean and Melchiorri 2002; 
Hannestad and Mortsell 2002; Melchiorri et al. 2003); CMB plus stellar ages (Jimenez et al. 
2003); and combinations of various databases (Kujat et al. 2002). Other potential methods for 
the determination of ojx have also widely discussed in literatures, such as the proposed SNAP 
satellite^ (Huterer and Turner 1999; Weller and Albrecht 2001; Weller and Albrecht 2002); ad- 
vanced gravitational wave detectors (Zhu, Fujimoto and Tatsumi 2001; Biesiada 2001); future 
SZ galaxy cluster surveys (Haiman, Mohr and Holder 2001); and gamma ray bursts (Choubey 
and King 2003; Takahashi et al. 2003). 

In this work, we shall consider the observational constraints on the DE equation of state 
parameterized by a redshift independent pressure-to-density ratio ojx arising from the latest ob- 
servations of angular size of high-z milliarcsecond compact radio sources compiled by Gurvits, 
Kellermann & Frey (1999) and the X-ray gas mass fraction data of clusters of galaxies published 
by Allen et al. (2002, 2003). The basics of a constant a>x assumption are two folds: on the one 
hand, the angular diameter distance used in this work is not sensitive to variations of a>x with 
redshift because it depends on ojx through multiple integrals (Maor et al. 2001 ; Maor et al. 2002; 
Wasserman 2002); on the other hand, for a wide class of quintessence models (particularly, those 
with tracking solutions), both of Q j and cjx vary very slowly (Zlatev et al. 1999; Steinhardt et al. 
1999; Efstathiou 1999), and an effective equation of state, Weff ~ / ajx(z)nxiz)dz/ / nx(z)dz IS a 
good approximation for analysis (Wang et al. 2000). We relaxe the usual constraint a>x > -1, be- 
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cause recent years there have been several models which predict a DE component with < -1 
(Parker and Raval 1999; Schulz and White 2001; Caldwell 2002; Maor et al. 2002; Frampton 
2003) and also we hope to explore its effects on the Wi determination. The confidence region 
on the (wj;, 0„,) plane obtained through a combined analysis of the two databases suggests 
-2.22 < (jji < 0.62 at 95.4% confidence level, which goes to -\ < oJx < 0.60 when the 
constraints,: > -1 is imposed. 

The plan of the paper is as follows. In the next section, we provide the bounds on from 
the angular size-redshift data. Constraints from the X-ray gas mass fraction of galaxy clusters 
are discussed in section 3. Finally, we present a combined analysis, our concluding remarks and 
discussion in section 4. Throughout of the paper, we assume a flat universe which is suggested 
by the measurements of the anisotropy of CMB and favoured by inflation scenario. 

2. Constraints from the angular size-redshift data 

We begin by evaluating the angular diameter distance 7>* as a function of redshift z. The redshift 
dependent Hubble parameter can be written as H(j) = HqE{z), where Hq = 100/j kms"^Mpc"^ is 
the Hubble constant at the present time. For a flat universe that contains (baryonic and cold dark) 
matter and dark energy with a constant <x>x (we ignore the radiation components in the universe 
that are not important for the cosmological tests considered in this work), we get (Turner and 
White 1997; Chiba et al. 1997; Zhu 1998) 



We first analyze the angular size-redshit data for milliarcsecond radio sources recently com- 
piled by Gurvits, Kellermann and Frey (1999) to constrain w, . The basics of the angular size- 
redshit test in the context of dark energy was first discussed in a theoretical viewpoint by Lima 
and Alcaniz (2000b) without using any database. They also provide an analytical closed form 
which determines how the redshift z,„, at which the angular size takes its minimal value, depends 
on Wj . Later on, using the same database compiled by Gurvits, Kellermann and Frey (1999), 
Lima and Alcaniz (2002) obtained 0„, ~ 0.2 and ~ -1.. A distinguishing characteristic of 
our analysis is that the usual constraints lo^ > -\ is relaxed. This database shown in Figure 1 is 
145 sources distributed into twelve redshift bins with about the same number of sources per bin. 
The lowest and highest redshift bins are centered at redshifts z = 0.52 and z = 3.6 respectively. 
We determine the model parameters and Q.„ through ax^ minimization method. The range of 
ajj spans the interval [-3,0] in steps of 0.01, while the range of il^ spans the interval [0, 1] also 
in steps of 0.01. 



where Oizf, ^m, (^x) = l/D^ is the angle subtended by an object of proper length / transverse to 
the line of sight and 6oi is the observed values of the angular size with errors cr, of the ith bin in 
the sample. The summation is over all 12 observational data points. 




E\z;n„,cox) = 0„(l+z)3+(l-Q„)(l+z)3(i^"^).(i) 




(2) 
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Fig. 1. Diagram of angular size vs redshift data for 145 compact radio sources (binned into 
12 bins) of Gurvits, KeUermann and Frey (1999). We assume the charateristic linear size 
/ = 22.64/j"'pc for theoretical curves. The soUd curve corresponds to our best fit with W;c = -1.19 
and £lm = 0.23, while the dashed and dot-dashed curves correspond to a A-dominated universe 
and the standard cold dark matter (SCDM) model respectively. 

As pointed out by the authors of previous analyses on this database (Gurvits, KeUermann and 
Frey 1999; Vishwakarma 2001; Alcaniz 2002; Zhu and Fujimoto 2002; Jain, Dev and Alcaniz 
2003; Chen and Ratra 2003), when one use the angular size data to constrain the cosmological 
parameters, the results will be strongly dependent on the characteristic length /. Therefore, in- 
stead of assuming a specific value for /, we have worked on the interval / = I5h~^ - SOh'^pc. 
In order to make the analysis independent of the choice of the characteristic length /, we also 
minimize equation (2) for /, o)x and simultaneously, which gives / = 22.64/j"'pc, (o^ = -1.19 
and Q.„ = 0.23 as the best fit. Figure 2 displays the 68.3% and 95.4% confidence level contours 
in the o)x) plane using the lower shaded and the lower plus darker shaded areas respectively. 
It is clear from the figure hat o)x is poorly constrained from the angular size-redshift data alone, 
which only gives cjx < -0.32 at 95.4% confidence level. However, as we shall see in Sec.4, 
when we combine this test with the X-ray gas mass fraction test, we could get fairly stringent 
constraints on both (o^ and 

3. Constraints from the galaxy clusters X-ray data 

Clusters of galaxies are the largest viriaUzed systems in the universe, and their masses can be 
estimated by X-ray and optical observations, as well as gravitational lensing measurements. A 
comparison of the gas mass fraction, /gas = Mgas/Mtot, as inferred from X-ray observations of 
clusters of galaxies to the cosmic baryon fraction can provide a direct constraint on the density 
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Fig. 2. Confidence region plot of the best fit to the database of the angular size-redshift data 
compiled by Gurvits, Keller and Frey (1999) - see the text for a detailed description of the 
method. The 68% and 95% confidence levels in the (Q™, a>x) plane are shown in lower shaded 
and lower + darker shaded areas respectively. 

parameter of the universe Qm (White et. al. 1993). Moreover, assuming the gas mass fraction is 
constant in cosmic time, Sasaki (1996) show that the /gas data of clusters of galaxies at difl'erent 
redshifts also provide an efficient way to constrain other cosmological parameters decribing the 
geometry of the universe. This is based on the fact that the measured /gas values for each cluster 
of galaxies depend on the assumed angular diameter distances to the sources as /gas «: [D^]^^^. 
The ture, underlying cosmology should be the one which make these measured /gas values to be 
invariant with redshift (Sasaki 1996; Allen at al. 2003). 

Using the Chandra observational data, Allen et al. (2002; 2003) have got the /gas profiles for 
the 10 relaxed clusters. Except for Abell 963, the /gas profiles of the other 9 clusters appear to 
have converged or be close to converging with a canonical radius r25oo, which is defined as the 
radius within which the mean mass density is 2500 times the critical density of the universe at 
the redshift of the cluster (Allen et al. 2002, 2003). The gas mass fraction values of these nine 
clusters at risoo (or at the outermost radii studied for PKS0745-191 and Abell 478) were shown 
in Figure 5 of AUen et al. (2003). We will use this database to constrain the equation of state 
of the dark energy component, tOx. Our analysis of the present data is very similar to the one 
performed by Lima et al. (2003). However, in additional to including new data from Allen et al. 
(2003), we also take into account the bias between the baryon fractions in galaxy clusters and in 
the universe as a whole. Following Allen et al. (2002), we have the model function as 



bQ.b \ h ^CDM^^') 



3/2 



(1 +0.19^1/2) [0.5 D^{zi\ux,0.,„) 



(3) 
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Fig. 3. Confidence region plot of the best fit to the /gas of 9 clusters published by Allen et al. 
(2002,2003) - see the text for a detailed description of the method. The 68% and 95% confi- 
dence levels in the ojx-^m plane are shown in lower shaded and lower + darker shaded areas 
respectively. 

where the bias factor b =^ 0.93 (Bialek et al. 2001; Allen et al. 2003) is a parameter motivated by 
gas dynamical simulations, which suggest that the baryon fraction in clusters is shghtly depressed 
with respect to the Universe as a whole (Cen and Ostriker 1994; Eke, Navarro and Frenk 1998; 
Frenk et al. 1999; Bialek et al. 2001). The term represents the change in the Hubble 

parameter from the defaut value of Hq = 50km s~^ Mpc"^ and the ratio D^(.j^^{zi)/D^izi; (o^, S^m) 
accounts for the deviations of the model considering from the default standard cold dark matter 
(SCDM) cosmology. 

Again, we determine cOx and through ax^ minimization method with the same parameter 
ranges and steps as last section. We constrain Clmh^ = 0.0205 ± 0.0018, the bound from the 
primodial nucleosynthesis (O'Meara et al. 2001), and h = 0.72 ± 0.08, the final result from the 
Hubble Key Project by Freedman et al. (2001). The;^'^ difl'erence between the model function 
and SCDM data is then (Allen et al. 2003) 

, /gas"^(Z(> ^i^a;) ^^m) ~ /gas.oi 
X\0jx, ^m) = 2_j 2 ■ 

/=1 /gas,/ 

where f^^izc, a>x, ^m) refers to equation (3), /gas.oi is the measured /gas with the defaut SCDM 
cosmology, and o-f^^^. is the symmetric root-mean- square errors {i refers to the fth data point, with 
totally 9 data). The sunmiation is over all of the observational data points. 

Figure 3 displays the 68.3% and 95.4% confidence level contours in the (ojx, ^m) plane of 
our analysis using the lower shaded and the lower plus darker shaded areas respectively. The best 
fit happans at (Ox = -0.86 and - 0.30 . As shown in the figure, although the X-ray gas mass 




Qfo/?^ - 0.0205 
0.0018 



h-0.12 
0.08 



(4) 
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Fig. 4. Confidence region plot of the best fit from a combined analysis for the angular size- 
redshift data (Gurvits et al. 1999) and the X-ray gas mass fractions of 9 clusters (Allen et al. 
2002, 2003). The 68% and 95.4% confidence levels in the ojx-^m plane are shown in lower 
shaded and lower + darker shaded areas respectively. The best fit happans at a>x = -1.16 and 
Q„ = 0.29. 

fraction data constrains the density parameter very stringently, it still poorly hmits the dark 
energy equation of state (Ox- The situation can be dramatically improved when the two databases 
are combined to analysis, in particularly, a nontrivial lower bound on cj^ will be obtained (see 
below). 

4. Combined analysis, conclusion and discussion 

Now we present our combined analysis of the constraints from the angular size-redshift data 
and the X-ray gas mass fraction of galaxy clusters and summarize our results. In Figure 4, we 
display the 68.3% and 95.4% confidence level contours in the (w^, Q^) plane using the lower 
shaded and the lower plus darker shaded areas respectively. The best fit happans at = -1.16 
and Q.m = 0.29. As it shown, fairly stringent bounds on both atx and Q.m are obtained, with 
-2.22 < ojx < -0.62 and 0.28 < Sim < 0.32 at the 95.4% confidence level. The bound on (Ox 
goes to -1 <cjx< -0.60 when the constraint cjx^ -1 is imposed. 

Although precise determinations of oJx and its possible evolution with cosmic time are crucial 
for deciphering the mystery of DE, currently cjx has not been determined quite well even with an 
assumption of ojx being constant (Hannestad and Mortsell 2002; Spergel et al. 2003; Takahashi 
et al. 2003). It is worthy of determining a>x using a joint analysis. In this paper we have shown 
that stringent constraints on oJx can be obtained from the combination analysis of the angular 
size-redshift data and the X-ray mass fraction data of clusters, which is a complementary to 



8 Zhu, Z.-H., Fujimoto, M.-K. & He, X.-T.: Determining ai^ from @-z and /gas data 

other joint analyses. At this point we compare our resuhs with other recent determinations of 
from independent methods. For the usual quintessence model (i.e., the constraint a>x > -1 
is imposed), Garnavich et al. (1998) found oj^ < -0.55 using the SNela data from the High-z 
Supernova Search Team, while Lima and Alcaniz (2002) obtained W v < -0.50 using the anugular 
size-redshift data from Gurvits, Kellerman and Frey (1999) (95% confidence level). Our result 
of a>x < -0.60 is a little bit more stringent than theirs. However Bean and MelchioiTi (2002) 
found an even better constraint, Wi < -0.85, by analyzing SNela data and measurements of LSS 
and the positions of the acoustic peaks in the CMB spectrum. For the more general dark energy 
model including either normal XCDM, as well as the extended or phantom energy (i.e., the 
constraint oj^ > -1 is relaxed), Hannestad and Mortsell (2002) combined CMB, LSS and SNela 
data for analyzing and obtained -2.68 < < -0.78 at 95.4% confidence level, whose lower and 
upper bounds are a little bit lower than ours (-2.22 < a>x < -0.62 at 95.4% confidence level). 
Recently, Schuecker et al. (2003) combined REFLEX X-ray clusters and SNela data to obtain 
-1.30 < a>x < -0.65 with Icr statistical significance. From Figure 4, it is found our Icr result 
is -1.72 < a>x < -0.83, which is comparable with the results of Schuecker et al. (2003). Using 
the X-ray gas mass fraction of 6 galaxy clusters, Lima et al. (2003) found -2.08 < a>x < -0.60 
(Icr level), which is less stringent than the result presented in this work. This is because we used 
more X-ray gas mass fraction data of galaxy clusters and combined the angular size-redshift data 
of compact radio sources for analyzing. The analysis presented here reinforces the interest in 
precise measurements of angular size of distant compact radio sources and statistical studies of 
the intrinsic length distribution of the sources. It is also hopefully that our constraints will be 
dramatically improved after more acurate X-ray data from Chandra and XMM-Newton become 
available near future. 
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